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ABSTRACT 

We apply a derivative expansion to the Legendre effective action ffow equations 
of 0{N) symmetric scalar field theory, making no other approximation. We 
calculate the critical exponents r], u, and uj at the both the leading and second 
order of the expansion, associated to the three dimensional Wilson-Fisher fixed 
points, at various values of N. In addition, we show how the derivative expansion 
reproduces exactly known results, at special values N = oo, —2, —4, ■ • 



1 Introduction 



In contrast to other analytic approaches, many different non-perturbative approx- 
imations can be envisaged within the Wilsonian renormahzation group frame- 
work |l| which preserve the structure of the continuum hmit {viz. renormahsabil- 
ity) 0. This fact alone makes a Wilsonian renormahzation group calculational 
framework potentially extremely powerful, and is surely in large part responsible 
for the resurgence of interest in these methods over the last few years. This frame- 
work holds other advantages also, for example enabling systematic searches for all 
(non-perturbative) continuum limits with the given field content and symmetries, 
and direct calculation of continuum quantities 0-0 {i.e. without introducing cor- 
responding bare quantities). Against this background, it is surely important to 
test and understand as thoroughly as possible the reliability and accuracy of the 
different possible systematic non-perturbative approximation schemes. 

Two such schemes appear particularly natural for the Wilsonian effective ac- 
tion: truncations in the operator content to a finite set of operators; or -less 
severely- an expansion and truncation in powers of space-time derivatives where 
no other approximation is made {a.k.a. "derivative expansion"). While carefully 
chosen truncations to a finite set of operators can be very successful in particular 
cases !§, 0! §1) there are reasons to expect these to have limited reliability and 
accuracy in general non-perturbative situations [0. It is also difficult to see how 
to keep such truncations systematic in practice beyond the lowest orders of the 
powers of derivatives they incorporate. In contrast, there are several reasons why 
one should expect derivative expansions to be well behaved @] at least in many 
situations of interest, while derivative expansions automatically yield a sequence 
of systematic approximations (since there is expansion in only one 'parameter'). 
Additionally, since a derivative expansion in effect incorporates at each new level 
of approximation a new infinite set of operators^, these expansions provide the 
limits of reliability and accuracy that can anyway be reasonably hoped for from 
truncations to any finite set of operators. 

The lowest order in the derivative expansion, ^0{d^)\ corresponds to the al- 
ready well established Local Potential Approximation in which the effective 
interactions are restricted to that of a general potential, and has proved reliable 
and reasonably accurate in a large variety of circumstances (see the references 
collected in 0, 0], and refs.^, ||],[0-[|16]). 0(9^) - second order in the deriva- 



tive expansion, has been explored in ||^,[|T3-[^. The derivative expansion for 
the Legendre effective action with infrared cutoff (this being the one-particle ir- 
reducible part of the Wilson/Polchinski effective actions ||2|) preserves more of 

^in bosonic theories, corresponding to all powers of the fields 
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the structure P, H {e.g. preservation of quantization of anomalous scaling di- 



mensions after derivative expansion, for certain cutoffs ||T9|, [Ty, as we review 
later). Very encouraging accuracy has been found in the results with this method 
for the non-perturbative massless quantum field theories of a single scalar field, 
corresponding to the two and three dimensional Ising model fixed points |jl9|, ^ , 



and the infinite sequence of two dimensional multicritical fixed points [2C]. These 
latter results are particularly significant since these calculations lie well outside 
the practical capabilities of all other standard approximation methods (includ- 
ing lattice methods) ||2^. Excellent results have also been found for the massive 
quantum field theory defined about three dimensional Ising model fixed point . 

Here we extend these second order calculations to the non-perturbative mass- 
less 0{N) scalar field theories in three dimensions, corresponding to the Wilson- 
Fisher fixed points We remind the reader that the following values of are 
experimentally realised in critical phenomena: = corresponding to polymers; 

= 1 corresponding to 'classical' critical liquid-vapour phase transitions, alloy 
order-disorder transitions, and uniaxial (Ising) ferromagnets; N = 2 correspond- 
ing to the He^ superfiuid phase transition, and to planar ferromagnets; = 3 
corresponding to (Heisenberg) ferromagnets; and A^ = 4 corresponding to the 



chiral phase transition for two quark flavours [^, ^ 



Our results for these cases, compared to the worlds best present estimates, are 
fair for the lowest order, and already rather goodQ at second order of approxima- 
tion. As A^ increases beyond four however, the results for t] particularly become 
poorer but we find good results for uj, until as A^ — >■ oo, we achieve the exact re- 



sults for all exponents 22]. We show also that at A^ = —2, —4, ■ ■ -, the results 



from derivative expansion coincide with the known exact results [0, |2g, |26 



At A^ = oo, we give explicit analytic solutions for the fixed point effective po- 
tential, and the eigenperturbations in the potential and kinetic terms. We obtain 
analytically all the corresponding eigenvalues, even though the perturbations to 
the kinetic term actually lie outside the exactly soluble sector in large A^ meth- 



ods [16 . At the special values A^ = —2, —4, ■ ■ -, we obtain analytically the known 



exact eigenvalues but also find intruiging hints for some novel fixed points. 

Most importantly in our opinion, we continue to find (for all A^), that the 
derivative expansion is qualitatively reliable - in the sense that no spurious fixed 
points or eigenvalues are found (in contrast to the generic situation for finite 
numbers of operators 0, |27|]). 

The organisation of the paper is as follows. In section 2 we review the method 



as set up in and provide the slight extensions necessary for the A^ compo- 



^with the exception of ui when > 1 
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nent case. In section 3 we develop the 0(9°) approximation. We remind the 
reader why the only approximation that actually need be made is to truncate 
the derivative expansion of the full connected two-point function. We derive the 
boundary conditions required for the fixed point and eigenperturbations. In sec- 
tion 4 we develop the 0(9^) approximation, again develop the relevant boundary 
conditions, identify the one redundant operator expected, and discuss the special 
behaviour at = 1. Section 5 deals with the case N = oo. We here derive the 
exact explicit solutions already mentioned above, and identify in passing, the line 
of Gaussian fixed points that appear at = oo, and the end-point responsible 
for the Bardeen-Mosher-Bander phenomenon. We point out here the possibil- 
ity to search for fixed points lieing outside the exactly soluble sector in large N 
methods, and provide physical justification for our analytic expressions for the 
eigenvalues for the kinetic term perturbations, which also lie outside the exactly 
soluble sector. In section 6, we show how to recover at = —2 and N = —4, 
the known exact results for certain eigenvalues. At the same time we find at 
A^ = —4, that there might exist a non-trivial fixed point with two exactly soluble 
irrational eigenvalues. A similar pattern is expected for the other negative even 
A^. Finally, in section 7 we present and discuss the numerical results. In deriving 
these, a number of numerical problems had to be understood and overcome, as 
described in appendix A. Appendix B provides some necessary beyond leading 
order asymptotic expressions. 



2 Review 

In this section we review, extending slightly where necessary, the derivative ex- 
pansion method developed in W^. We take the following as our partition function. 



expW[J] = y'l)v^exp|-iv^.C-\v^-5Ao[v?] + J.<^} ■ (1) 

We use a condensed notation, so that two-point functions are regarded as matri- 
ces in position or momentum space, as well as any internal index space; one-point 
functions as vectors; and contractions indicated by a dot. We will work in D Eu- 
clidean dimensions with a real scalar field with A^ components: (f"", a = 1, ■ ■ ■ , N. 
Sao is the full bare action for the theory, and is taken to have an internal 0{N) 
symmetry. C = C(g,A) is taken to be an additive cutoff for convenience jTol] , 
satisfying C — * as g — ^ 0, and q'^C{q, A) ^ oo as g ^ oo. 

From (|l]) we have 
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Re- writing this in terms of the Legendre effective action, T\, 
T/^[ip'^] + ^ip'^ .C~^ .ip'^ = —W[J] + J.f'^, where the classical field ip'^ 
gives 



d_ 

dA 



1 



-tr 




1 + C. 



defined by 
= SW/SJ, 



(3) 



The trace represents a sum over the internal spin indices and an integration over 
momentum/position space. As in ref. it is convenient to write the trace as an 
integral over momentum space and factor out the D-dimensional solid angle. At 
the same time, we drop the superscript on the classical field: 



dA 



-tr / dq 
2 Jo 



q^~^ dC{q,A) 
C{q,A) OA 



1+C. 



(4) 



where Q = 2/[r(D/2)(47r)^/2] jg 

the solid angle of & D — 1-dimensional sphere 
divided by {27i)^, the brackets (■ ■ ■) represent an average over all directions of 
the momentum q, and the trace now represents only the trace over the internal 
indices. 

At a fixed point, the field (p scales anomalously as ~ A'^'^, where = 
^{D — 2 + 7]) and t] is the anomalous scaling dimension. Considering the defining 
expression for Fa, ^a['p] = —\'p.C~^.'p — Wa[J] + J.'P, dimensional analysis then 
shows that we require C to behave as follows. 



C(g,A) ^ A''-2c'(g2/A2) 



(5) 



for some C, if we are to express Fa as independent of A at the fixed point. From 
now on we write C(g, A) as A^~^C*(g^/A^) and drop the tilde on the scaled C. As 
we are only interested in the behaviour near or at fixed points this is a sensible 
definition and we can take to be a constant. This is the reason for choosing 
an additive cutoff: a multiplicative cutoff Cjr cannot scale homogeneously (to 
absorb A''), and remain consistent with its normalisation in the ultraviolet {viz. 
Cm — > 1 as q — > oo ) 111. Now it is necessary (for Fa to be independent of A at 
fixed points) to re-write the equations in terms of dimensionless variables, 

q Aq 

(^(Aq) ^ A^-^-y,(q) , (6) 

and we define t = log(/i/A) (where /i is an arbitrary finite energy scale). Finally 
we rescale the fields and the effective action to absorb the factor of f2/2 in (^, 
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where C is a numerical normalization factor to be chosen for later convenience. 
Upon doing this we get, 



(7) 



Ctr / dqq 



D-l 



' q dCjq^ 
C{q^) dq 



+ 2-7] 



1 + a 



-1 



iq, -qj 



In the above — '^■j^ is the field counting operator: it counts the number 
of occurrences of the field in a given vertex, and arises due to the scaling of 
the field in (^. Ag is the momentum counting operator plus the dimension of 
space D and arises through the rescaling of the momenta in equation (^. It can 
be represented as 



Operating on a given vertex it counts the total number of derivatives acting on 
the fields ip. Equation (|^) will be the starting point for all the work from now 
onwards. Notice that for the first time t] explicitly appears in the equation. 

We write Tt[ip] as the space-time integral of an effective Lagrangian expanded 
in powers of derivatives, 

= Jd^'x {v{v\ t) + \k{^\ t) {d.^'^f + \z{^\ t) {^^d.v^f + ■ ■ j (9) 

Each linearly independent (under integration by parts) combination of differ- 
entiated fields will carry its own general (t-dependent) coefficient. The global 
0{N) symmetry forces us to choose the coefficient functions to be functions of 
ip"^. We will require that the fixed point solutions for V,K,Z etc. will be non- 
singular for all <y9^, that perturbations about these solutions grow no faster than 
a power P^, ^, ^, and that K{0) ^ 0. It will be seen that the imposition of power 
law growth results in a quantized spectrum. 

It was realized some time ago that a general cutoff function C(g^), for example 
based upon an exponential function, leaves rj undetermined and dependent upon 
an unphysical parameter |T^, ^. However, dimensional analysis indicates that 
if C{q^) is chosen to be homogenous in q^ then equation (J^) is invariant under 
a global rescaling symmetry |TP|. This rescaling symmetry overdetermines the 



equations, so that solutions only exist for discrete values of rj. We will follow this 
path and take C{q^) = q"^^, for k a positive integer. From general arguments it 



can be shown that the convergence is slower the higher the value of k [Q, p!9[] . 
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On the other hand, if we wish the integrals in equation (0) to be ultra-violet 



convergent we must take k > D/2-1 [0. Therefore we take k to be the smallest 
integer greater than D/2 — 1, so for D = 3we take k = 1. The rescaling symmetry 
is made explicit by choosing the following (non-physical) scaling dimensions, as 
follows from (|^) and the definition of Fj. 

[q^'] = 1 [d''] = 1 [^^] = k + D/2 (10) 

[V] = D [K] = -2{k + 1) [Z] = -2{2k + 1) - D . 

The expansion is performed by substituting equation (^) into equation (0) and 
expanding the right hand side up to a maximum number of derivatives. The 
angular average in (|^) can be easily computed by translating them into invariant 
tensors, e.g. {q^q'^) = q^d^^ / D etc. From now on we will specialize to D = 3. 



3 Lowest Order of Approximation 

To lowest order in the expansion we drop all the derivatives from the right hand 
side of (|^. At the same time, a reasonable ansatz for Fj would be to keep 
only a general potential V{^p^,t) in (^, setting K = 1, and Z etc. to zero. 
It is worthwhile to remark however, that there is no need for a further ansatz 
beyond that of dropping all derivative terms from the right hand side of (|^), 
i.e. from the full connected two-point function (and similarly at higher orders 
of the derivative expansion): the form of Ff is then determined by consistency 
arguments []Ty|. Substituting (^ into (|^), we see that the coefficient functions 
K,Z, etc. are determined by linear equations given by the vanishing of the left 
hand side of (|^). Consider the equation for 

^ir(^^ t) + (1 + V)V^K\^\ t) + T^K{^\ t) = (11) 
where ' = -g^sj- We see that at fixed points, where K is independent of t, 

If K{ip'^) is to be non-singular (at = 0) and -ft'(O) 7^ 0, we must have 77 = 0. 
Therefore K{if'^) is a constant. Using the scaling symmetry we can set this 
constant to be iiT = 1. If we now consider the equation for Z we get 

^Z(y,2, t) + ^'Z'{^\ t) + Z{^\ t) = (13) 

Hence we see that at fixed points, Z will satisfy 

Z OC {^'Y' (14) 
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To obtain non-singular behaviour, we must set the constant of proportionahty to 
zero and thus = 0. Considering the form of a general term in the expansion 

we see that a similar conclusion must hold for all other terms. A general term, 
if, will satisfy 

j^H{^\t) + ^^H\^\t) + n^H{^\t)+naH{^\t)-?,H{^\t) = Q , (15) 

where n^p is the number of 9?'s occurring in the term multiplying if, divided by 
two and uq is the number of derivatives occurring in the term multiplying H . 
Hence we see that at a fixed point, we have 

H OC (y,2)-K+na-3) _ ^^g) 

Thus for terms of higher order in the expansion, where > 2 and ng > 2, H 
will be singular unless the constant of proportionality is set to zero. Therefore, 
at leading order, the derivative expansion must reproduce the local potential 
expansion at fixed points, and thus we take 

r,M = jd''xv{ip\t) + \{d,ip^f . (17) 



Since we discard all derivatives from 



OlpOip 



(q, — q) , we have that ip is 



effectively a constant and the operator diagonal in momentum space (as explained 
further in the next section), thus it is the matrix inverse of 



5"* + C{q' 



6ip''{-q)6ip''{q) 



The inverse and trace in (0), are now easily accomplished by noting that a matrix 
A6°''' + Bip'^ip'^ has — 1 eigenvalues A and one eigenvalue A + ip^B: 



dt 



+<^v(v9^t)-3l^(<^^^) 

N -I 



(19) 



Jo^^^ l + 2V'{ip^,t)q^ + q^ ' 1 + {2V'{ip'^,t) + Acp^V"{ip'^,t))q'^ + q"^ 



+ 



Performing the g-integrals, and setting ( to l/27r yields the equation at leading 
order, 



dt 



+ ip'V'{ip\t)~3V{v',t) 



1 



^2 + 2V'{ip^,t)+4:ip^V"{ip^,t) 
N-1 



2 + 2V'{^^,t) 



(20) 



Let us now discuss fixed point solutions of (POD, that is solutions with dV/dt = 
0. At first sight it may seem that (^0]) has infinitely many fixed point solutions. 
In fact this is not the case, as only finitely many solutions do not end in a 
singularity |]14|, Of course this is sensible on physical grounds, as the fixed 
points correspond to massless continuum limits {i.e. continuous phase transitions) 
with the prescribed field content. Considering the form of (pO|) we see that either 
V{ip'^) is only defined for yj^ < yj^ and ends at a singularity, at ip"^ = ifl as follows, 

(where we have suppressed non-singular and lower order singular parts), which 
we disregard as unphysical, or it is defined for all > 0, and in particular 
for if^ oo. In this regime, we find that either V{(.p^) is just a constant (the 
Gaussian fixed point), or V{'^'^) must satisfy the following, 

^„„.)3^(^,_L_i)_L_,o((.V) (22) 

for some constant Ay. If we consider linear perturbations about this solution we 
see that for large ip^ the perturbation must behave as a linear combination of 
(v?^)^ and exp(|[3074^]^/^(93^)^). The first perturbation merely alters the value of 
Ay, whilst we disallow the second as it is incompatible with (^). We see that, 
apart from the Gaussian fixed point, the solution space of fixed point solutions 
defined for all yj^, forms an isolated one parameter set, parametrized by the value 
of Ay. Since we require V{(.p^) also to be non-singular at the origin, setting = 
in equation ( pOD we obtain, 

3^(0) = -— ^== (23) 
12 + 2\/'(0) 



We now have a second order differential equation with two boundary conditions: 
([2^j23|) . Thus we expect at most a discrete set of acceptable solutions. In fact 
we only find two:- the Gaussian fixed point, V{lp^) = and an approximation 
to the Wilson-Fisher fixed point To find this non-trivial solution we need to 
rely on numerical methods, as outlined in the appendix A. The results of these 
calculations are shown in figure |I], for various values of A^. 

To calculate the critical exponents, we linearize about this fixed point poten- 
tial V = V*{ip^). Writing V{ip'^,t) = V* (if^) + 6V {ip"^ , t) , where 6V{ip^,t) is given 
by ee'^^v{ip'^), with £ <^ 1, and expanding to first order in e, we have 

(A-3)t;((^ )+^ v[^)= . ...3/2 + . . ■ ...... ...3/2 (24) 



(2 + 2V'{^^)f^ (2 + 2V'{^^) + Vv//(y,2)) 



8 




-2.0 ' • ' • ' • ' • ' 

0.0 1.0 2.0 3.0 4.0 

Figure 1: The potential at leading order in the derivative expansion. 

Requiring that f (v?^) is not singular at the origin ip'^ = 0, we obtain 

TV 

''-^^"^- (2 + 2^0))^/^ "''°^ ■ '''' 

Thus for 7^ we may choose, by linearity of the perturbation, the normalisation 
condition v{0) = 1. For = 0, (p5|) implies that either A = 3 or v{0) = 0. We 
discard the case of A = 3 as this corresponds to the uninteresting case of the 
vacuum energy operator f (v?^) = 1- Therefore for = one has v{0) = 0, and 
by linearity we can take v'{0) = 1. Thus for any iV > 0, we obtain two boundary 
conditions at the origin. For large we see that the perturbation will be a linear 
combination of 

i^y-' (26) 

and exp(|[30y4^]^''^(v3^)^). Once more we will enforce a coefficient of zero on the 
latter perturbation, as we require the perturbations to grow no faster than a power 
in y?^ [0, ^, ^. The imposition of this condition will ensure that the solutions form 
an isolated one parameter set. Upon also imposing the two boundary conditions 
at the origin, we expect at most a discrete number of solutions to the eigenvalue 
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problem for A, which is indeed the case. We find only one positive eigenvalue, 
which yields v through v = 1/A. The least negative eigenvalue yields the first 
correction to scaling through uj = —A. The results are summarized in the table 
presented in section 7. 



4 Second Order of Approximation 



The equations at the second order in the expansion are calculated by substituting 
expression into equation (|^ and dropping all terms with more than two 
derivatives from the right hand side. Following a similar argument to that given 
in the previous section, we see that terms with more than two derivatives are 
forced to vanish, hence at second order of the derivative expansion, only the 
terms explicitly written in (P) survive. A little bit of thought shows that a similar 
conclusion holds at all higher orders of the approximation: if we substitute an 
expression into the equation that is of higher order than the order we are working 



at, then we are forced to set the higher order terms to zero |T^ 



The computation of the inverse in equation (0), is not as straightforward as in 
the leading order case. We regard [1 + C.6'^Tt/6ip6ip]~^ as a differential operator: 



C. 



n -1 



where 



iq, -qj 



Q 



Jd^xQ 

1+C. 



6^^ + C. 



H -1 



ix,yje 



-iq.x 



n -1 



(27) 
(28) 



Q is a function of q and V'(x) and its derivatives evaluated at x, and C and 
g^^F^b are written as differential operators: C(— □), 



ab/ 



-2^\d,^')Zd, - 2^^^\d^^'fZ' - 5^\v^U^')Z - v\U^')Z 



2^\d^v'){^'d^^')Z' - 2^'^^\^-d,v-)Z'd, - 2<^V(<^'5m<^')'^" 



(29) 



Define z/ as the expression obtained by dropping all terms containing differentials 
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of Lf from C(g^)e ''^ 



2^„— jq.x ^■^r „iq.x 



ra6T // 



ab ) 



(Recall C(g^) = in the present case.) Noting that from 



we see that Q"^ satisfies the following expression ||19|| , 



1 + vU + {l + v) 



-1 I , ,cd/^db 



-iq.x 



C{-u) 



(30) 



(31) 



(32) 



(as may be verified by multiplying by + z/^"). By construction, the term in 
curly brackets in (|3^) is at least of first order in the derivative expansion. The 
derivative expansion to any order may thus be derived by iterating expression 



(p^), starting with the zeroth order Q 



ab 



[i + y) 



We iterated just to two derivatives [using the fact that 1 + i/ is of the form 



^ab ^ ^ab ^ ^^ab ^ Bip^ip\ gO (1 + Z/)-^^ = ^ 



B 



AiATW)^"^ ' identity 
— (1 + v)~^{d^v){l + ]. Even so, this is a long calculation. 



which we performed using the symbolic manipulation package Form ||2^ . Finally, 
in preparation for the N = oo limit, we scale the equations as follows p2|, p!6|. 



ip = VNif, V = NV, K = K, 
From the zeroth order part of Q"""", we obtain 



Z/N . 



^V + {l + v)v>'V'-3V 



(1-V4)x 



(33) 
(34) 



1 



+ 



1 - 1/N 



(where here and from now on we drop the tildes on yj, V", K and Z). Similarly, 
from the [df^ip'^)'^ part of Q"" we obtain, 

d_ 

dt 



K+{1 + 7])ip'K' + r]K 



TT 



-'4,2,1 I - — — 1- : 



3 TV 



+ -'4,1,2 I - — — 4 — r 



3 N 



N 



3 N 



16 h^2,iV" K' if ^ fl6V"K'ip 16 V" Zip' 



2,0,2 



K' 



'2,2,0 



~iV^^ TV 



(35) 
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where for compactness we have defined the integrals 



j,k,m 



[1 + {W + Vy")g2 + (is: + i{)'^Z)q^y (1 + Wq^ + Kq^) 



(36) 



These integrals can be evaluated in closed form, but (|35|) would then be too long 
to display. Furthermore, as discussed in appendix A, it proves numerically not 
sensible to do so. Even in this compact form however, the fiow equation for Z 
(which follows from the {ip'^d^ip'^Y part of Q"") is too long to display p9| , 

d 



dt 



Z + {l+r])ip^Z' +{1 + 27])Z 



(37) 



As in the previous section, we obtain the boundary conditions for fixed point 
solutions, from the requirement that these solutions be non-singular. Thus for 
large 99^, either the solution is the trivial Gaussian fixed point [rj = 0,V = 
l/3^/2, K = 1, Z = 0), or the solutions behave as follows, 

Viip^) ~ A, {ip^)TT^ + ... (38) 
Kiip^) ~ Ak {ip^)TT^ + ... (39) 
Z{ip^) ~ A, {ip^yTW + ... . (40) 

We can also study the perturbations about these solutions, and forcing that 
the perturbations grow no faster than a power, we see that solutions satisfy- 
form an isolated four-parameter set, that is including rj. From 
the requirement that the solutions are non-singular at the origin, substituting 
= in (^) - ([37|), we now obtain three conditions. Using the scaling sym- 
metry (|ll|), we can impose one further condition. We take K{0) = 1 (with other 
possible solutions being reached by using the reparametrization invariance). Thus 
with four conditions imposed on a four-parameter set, we expect only a discrete 
number of solutions, in particular for 77. In this way, (p^) - (pTj) at fixed points, 
form non-linear eigenvalue equations for t]. 

Again we only find two solutions:- the Gaussian point mentioned above and an 
approximation to the Wilson-Fisher fixed point. The results for 77 are summarized 
in the table presented in section 7, and the results for the fixed point solutions 
shown in figures ^ to ^. 

To calculate the other critical exponents (at second order of approximation) 
we write. 



mg 



to (40 
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0.0 



0.0 



0.5 



1.0 



1.5 



Figure 2: The potential at second order in the derivative expansion. 



(41) 
(42) 



Z{if\t) = Z*{ip^) + 5Z{ip\t) 

= Z*{^'')+ee^'z{^'') (43) 

where V*{ip'^)^ K*{(p'^), Z*{ip'^) (and rj*) are the fixed point solutions calculated 
above, and expand the equations to first order in e. If we insist that v,k and z 
grow no faster than a power at large ip'^ ^ Q , then asymptotic analysis shows 
that V, k and z grow according to their scaling dimension, 



z{^') 



--A 



l + T, 



--A 



1+2?) 



0^(V5 ) + 



(44) 
(45) 
(46) 
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N = 2 

N = 3 
N=4 
N = 10 

N = oo€ 



1.0 



-e- 



0.9 




0.8 I • ' • ' • 

0.0 0.5 1.0 1.5 

Figure 3: The K component at second order in the derivative expansion. 

and form an isolated four-parameter set (including A). As before, another three 
boundary conditions come from substituting (yj^ = 0, assuming non-singular lim- 
its. Using linearity to set v{0) = 1, we then expect at most a discrete set of 
solutions, which is what is found. Again we find just one positive eigenvalue, 
which yields u, and determine the least negative. These are shown in the table, 
presented in section 7. 

It is important to recognize that there are other solutions of the equations 
which do not correspond to critical indices These solutions are known as re- 
dundant perturbations and the eigenvalue corresponding to the solution depends 
on the exact form of the renormalization group chosen. The redundant pertur- 



bation reflects invariances of the equations [^, [T9[. In fact we expect and find 
only one redundant perturbation corresponding to the reparametrization invari- 
ance (pTf). This has eigenvalue A = 0, and (not normalised) 

V = -5y?V + 3V, k = -bif'^K' - AK, and z = -bif'^Z' - 9Z . (47) 

An obvious question is why no redundant perturbations are found at the lead- 
ing order. The answer is simple: the choice of i^' = 1 [viz. for all t) breaks the 
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Figure 4: The Z component at second order in the derivative expansion. 



reparametrization invariance, thus we should not expect to see any redundant per- 
turbation. Similarly we should not expect to see any perturbation corresponding 



to (^-translation invariance [pO| , pO| as the choice of V, K and Z being functions 
of breaks this invariance. 

The case of = 1 is special. Of course in this case there is only a discrete Z^ 
internal symmetry, but also the derivative expansion (Rl) now becomes, 



d^'x \V{^\ t) + i {k{^\ t) + ^^Z{^\ t)) {d,^f\ (48) 



We see that Z and K no longer have a separate significance, and we should 
consider only the function n = K + (p^Z. In fact at = 1, the Z component of 
the fixed point solution to (P^)-(P7|) is singular, diverging as ~ l/y?^ as ^ 

I5D+(^2(|3^), we find that at = 1 



(c/. appendix A). Taking ([3^), and forming 
all appearances of K and Z can be absorbed into /t, after which they coincide 



with the already analysed equations in ref. [IS], as expected. 
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5 Exact Results at = 00 



It is known that in the large hmit, the local potential approximation can 



effectively become exact |TT|, ^, [T^, |T6[, in the sense that {\L7\) then coincides with 
an exactly soluble subsector of the full flow equations. In this section we will 
show that the derivative expansion reproduces this result, and give an explicit 
expression for the non-trivial (Wilson-Fisher) fixed point V*{(f'^) and for the 
perturbations about this fixed point. These solutions play an important role 
for us, because we obtain from them numerically, by relaxation (c/. appendix A), 
the solutions at finite N. 

Taking the N = oo limit of equations ( |51D - (0), we obtain, 

|y + (l + ,)v.V-3V = - (49) 

Vk {2V' + 2VKj 

^K+il + ,)^^K' + ,K = (50) 

^/K [2V' + 2VK) 

^Z+(l + 77)<^2^' + (l + 2r/)Z = ■■■ . (51) 

The equations for K and V have decoupled from the equation for Z, and all the 
equations are now first order equations. 

The exactly soluble large N limit arises when the interactions in Ft (or the 
bare action S'aq) are functionals only of ip"^ |]16[. Evidently, this is true of V in 



(H), while all such derivative squared interactions can be cast (by integration by 
parts if necessary) in the form | {d^ip^)^ ^i'-p'^i i)- This coincides with the Z term 
in (^ but rules out a non-trivial K. Therefore we require^ K{(p'^,t) = 1. This is 
indeed a solution of (^OD, but only if we also set rj = 0. 



Note that if we take rj = 0, then from (|50| ) the fixed point solution has to be 
trivial: K{ip'^,t) = 1. It would be very interesting to perform a thorough search 
of the two parameter space {e.g. parametrized by V{0) and rj) of t-independent 



solutions to ( ^pO[ ), looking for non-singular non-trivial solutions, which would 
if they exist, thus correspond to novel fixed points with non-trivial K and rj 0. 
From the above analysis, such solutions for K would lie outside the exactly soluble 
sector of the large limit [0. 

We will see that non-trivial perturbations SK around the Wilson-Fisher fixed 
point, which (from the above) also lie outside the exactly soluble sector, yield 
very sensible results for the eigen-value spectrum. 

3 



We are free to choose the normalization to be unity by the scaling symmetry (hi 
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Setting r] = 0, and Ki^ip^) = 1, we have that the fixed point potentials satisfy, 



(^V -3V = ^ (52) 

Differentiating with respect to y?^ and setting Wlip"^) = V^'(y9^) we then have, 

W 

This is trivially solvable as a differential equation for yj^ with respect to W, 

ip^ = ±AVw + ^ViTw(i + -^) (54) 

where A > is a constant of integration. A > 1/ \pl corresponds to the line 
of Gaussian fixed points that appear at = cxd [|15], |16|, with the negative sign 
root for Lp^ < l/2\/2, and the positive sign root for ip"^ > 1/2a/2. A = 1/V2 is 
the end-point where the Bardeen-Mosher-Bander phenomenon appears ||3l|, |15[] . 
< A < I/V2 is not allowed, because there exists a range of field values (p"^ < ipl 
for which no real solution W exists. The value A = is the only remaining 
possibility, and is the solution we take since this corresponds to the Wilson- 



Fisher fixed point [11, 15, 16|. Inverting (B3) then gives two solutions, only one 



of which satisfies the original equation (Fl) however: 



W = {ip^f -l/2 + ^^^l + {^^Y . (55) 
Integrating this solution and comparing with (^), gives the potential: 

V = -\^' + \{^'f + \{{^'f + lY'" . (56) 

Substituting Ki^ip"^) = 1 and r] = into (|5TD, we find 
2v/ , V _ , 1024 ^ ,..4,^o^2 



TT L 6 

32 

+ — ip^ (4 /6,2,3 - 8 h,2,3 + ^4,2,2 + 4 /2,2,3 - 5 /2,2,2 + /2,2,l) V"V 

64 

+ _^2 (^^Q j^^^^^ h,2,3 - 11 /3,2,2 + 8 /3,2,3 - 8 /3,3,2 + 8 /3,3,1 + 8 15,3,2) V 

/256 256 256 416 2 ^ , 

+ I — ^ Z i6,3,2 ^ V Z 74,3,2 + V ^ -'4,2,3 ^ Z 74,2,2 

256 2n2 72 r 128 256 128 64 

H — ^ IV^ j ^ -'6,3,2 H — ^ -'6,2,3 ^ -'4,2,3 ^ 64 i4,2,2 H — ^ -'2,2,3 — ^ -'2,2,2 
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//3 



64 128 64 128 128 256 80 

^ -'2,2,1 H ^ -16,0,5 + -'6,3,2 H ^ -14,3,1 ^ -'6,1,4 H ^ -'4,1,4 + -'4,1,3 

128 16 64 64 128 112 

H ^ -'2,0,5 + -'2,0,3 ~" 32 /2,o,4 + -'2,3,0 + -'2,3,2 ^ -12,1,4 H ^ -12,1,3 

16 128 32 256 128 256 2 ^ . 

+ ~^ -'2,1,2 — -12,3,1 ^ -'4,0,4 — -14,0,5 ^ -'4,3,2 H ^ ^ -14,3,1 

256 2 \ 2 2y,r 7T 40 2 7/ , 64 2 

Z /6,2,3 ) V/ + ( - V5 Z /3,2,1 + 3 ^ ^5,2,2 - y Z /3,2,2 - y (/^ Z J5,2,3 

7T 7T ^^'^ 7T ^^7T ^ ^2 2 7, , ^ . ^ , 
^ Z J3,2,2 + 2 ^ -'3,2,1 ~2~ -'3,2,3 ^ ^ is, 2,3 ~'~ "3" "'7,2,3 + o Z i3,o,3 

32 „ . 32 



+ y Z /7,2,3 + {^'Z' /3,2,3 + ^ ^' ^5,2,2) } 



where the Ij,k,m were defined in (^61). Needless to say, the above repesents a great 
simphfication compared to the full Z equation [^]. An analytic solution of this 
equation might be possible, by the methods of ref.[|l^], but the numerical solution 
is straightforward. Imposing the constraint that Z exists for all ip^ we get the 
solution shown in figure ^ 

To calculate the critical exponents we again linearize about the fixed point by 
writing V{Lp'^,t) = V*{if'^) + ee^'^v{(f'^) and K{if'^,t) = K*{(f'^) + ee^^kiif"^), where 
V* and K* denote the fixed point solutions found above. From (|^) and (|50D, 

X ^ 2/ ^ I k + 2v' ,1 k 

^- + ^--^- = 2 (2 + 21^)3/2 + 2 7f + W ^'^^ 

- (^tI^ ' ^''^ 

where W is given in (^). These equations are straightforwardly soluble: e.g. 
from ( |53D and (0) we have by inspection^ 

k = a,{W/2)-'/' , (59) 

while substituting v = k{ip'^)f{ip'^), one finds 

(These expressions may be further simplified by defining z = e*^™*^ then 
W = — 1 and the multiplier of /c/4 above, becomes — 3z + 1/^.) 

Noting the simple zero of W{ip'^) at (f"^ = l/2y/2, we see that (|59D and (j^) 
are non-singular if and only if: either ak = k = and A = 3, 1, —1, ■ ■ ■, the exactly 

^following an observation by Haruhiko Terao 
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soluble potential perturbations and spectrum found previously [0, ||, |T^ ; or 
ciy = and A = 0, —2, —4, ■ ■ -. The A = solution of this set is the redundant 
perturbation (|i7|). The others lie outside the large N exactly soluble sector (as 
discussed above) and are thus presumably only approximate, however we see 
that their spectrum can be understood, as with the potential perturbations, by 
assigning scaling dimension [y?^] = 2 to y?^, regarding this as a composite field. 
This interpretation seems reasonable in view of the explicit appearance of such a 
field with this dimension in traditional large approaches p2 . 



6 Exact Results at iV = -2 



The fact that = —2 is also an exactly solvable case was noted long ago ||2^, ^ . 
The initial hope was to be able to combine the known results at = oo with 
those at = —2, by using Pade approximants, to gain information about the 



physically interesting cases A^ = 0,1,2,3 p^, ^ It was found that when 
A^ = —2, the exponents are always Gaussian, rj = and u = 1/2. We will 
show how the derivative expansion reproduces these results at the leading order 
of the approximation. At the end of the section, we also outline some interesting 
behaviour we observe at A^ = —4. 

Differentiating (^) with respect to ip^, and setting ip^ = we have that fixed 
point potentials V^ip"^) satisfy, 

- 2V'(0) = ^^A../o ^"(0) (61) 

^ ' {2 + 2y'(0)}3/2 ^ ' ^ ' 

Thus at A^ = -2, V'{0) = 0. Substituting this into (|0D, we find 1/(0) = -^2/3. 
Therefore in this case, V{0) and V'{0) are exactly determined, but V"{0) is an a 
priori free parameter. Numerically, we find that V"{0) is fixed by the requirement 
that V{ip^) be non-singular for all ip"^ > 0. We find two solutions: the trivial 
tricritical Gaussian fixed point V{(p'^) = —y/2/3, and a non-trivial Wilson-Fisher 
fixed point with V"{0) > 0, satisfying the required asymptotic conditions at large 
ip"^ {cf. equation (p^)). 



Now consider the eigenperturbations. Differentiating ( p4D with respect to ip , 
and setting 99^ = we have 

(\ 9W(n^ + 2)^"(0) JiV + 2)t;-(0)V-(0) 

(A-2).(0)=^2V'(0) + 2}3/2-^ {2n0) + 2p/2 • 

Hence when A^ = —2, A = 2 or v'{0) = 0. For the case A = 2, we checked 
numerically that there exists a non-singular solution with the required asymptotic 
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behaviour This eigenvalue yields v = 1/2. For the case f'(0) = 0, setting 

= in (p4D we find (A — 3)f (0) = 0. Thus either A = 3, which corresponds to 
the trivial vacuum energy perturbation f (v?^) = 1, or f (0) = 0. 

In the latter case we now have f (0) = v'{0) = 0, so for non-trivial solutions of 
this type we may normalise v"{0) = 1. We expect from Sturm-Liouville theory 
that about each fixed point, for an infinite discrete set of eigenvalues A < 2, 
these solutions have the required asymptotic behaviour p^).p| The perturbations 
about the Gaussian fixed point can be determined exactly: it is easy to check 
from (0), that the eigenperturbations with the required asymptotic behaviour, 
correspond to the Laguerre polynomials [|l^, |15[ v{lp'^) = 2L^/^~^(y9^v^) (given 
by the formula L^{y) = ^eS'?/-"^ (e-^?/"+") fS^) with A = 3 - n, n = 0, 1, • • •. 
For n > 2, these indeed satisfy f (0) = f '(0) = and v"{0) = 1. 



As pointed out in Fisher |2^ we should also expect to find exact solutions 
for = —4, —6, ■ ■ ■. For = —4 we find, in a similar way to above, that 
V"'{0) is now an a priori free parameter, while either V"{0) = V'{0) = and 



V{0) = -23/2/3, or V"{0) = 12^/^5^/^ V'{0) = 1/5 and ^(0) = -^20/27. We 
have not checked however, that non-singular solutions actually exist with these 
boundary conditions (apart from of course the Gaussian solution V"'{0) = 
allowed by the first set). The first set of boundary conditions correspond to 
Fisher's findings since we then obtain (A — 2) (A — l)f'(O) = 0, so that the 



Gaussian-type eigenvalues A = 1, 2 are recovered when acceptable solutions exist 
with v'{0) ^ 0. (Additionally it may be shown that v"{0) = if A 7^ 1 or 2.) 
The second set leads to (6A^ — 9A — 10)i;'(0) = 0, suggesting the existence of a 
non-trivial fixed point with two exactly soluble eigenvalues A = 3/4 it -\/321/12. 

To confirm this finding (at this level of approximation) would require how- 
ever to check that a non-singular fixed point solution actually exists, as men- 
tioned above, and that eigenperturbations satisfying ( p6|) for large ip"^, exist 
with the above eigenvalues. Equally, we expect to recover Fisher's results at 
A^ = —6, —8, ■ ■ ■, and further novel candidates for fixed points. We do not pursue 
these interesting questions further. 



7 Numerical Results 

In the table, we display the numbers only to the accuracy necessary for compari- 
son to the worlds best estimates. We did not find second order results for A^ = 0, 

^This is not at variance with |^ since the source term J -tp in (|l]) generates 'anisotropic' |^ 
correlation functions. 
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.93 


.89 


100 


.00034 
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.994 


.998 


.989 


.991 
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.98 



Table 1: Critical exponents of the three-dimensional Wilson-Fisher fixed point. 
The first two orders of the derivative expansion are compared to a combination 
of the worlds best estimates with their errors, where available. 0{d^) is 

the leading order in the derivative expansion, where rj is identically zero for all 
A^, and 0(9^) is the second order in derivative expansion. 



for numerical reasons, as explained in appendix A. 

Since, as we have seen, already the Local Potential Approximation yields the 
known exact results at = —2, —4, ■ ■ -, and at = — oo, at the very least these 
approximations provide physically motivated interpolations between these exact 



values, which thus go beyond the early hopes of using Pade approximants p4 



In fact there is quite an impressive agreement between 0(9°) and estimates 
by other methods. Compared to these other methods, at first there is a gradual 
decrease in the accuracy of the prediction for z/, as N increases, whilst the pre- 
diction for UJ shows a slight improvement, until as — > oo all exponents become 
exactly determined. The results compare favourably with those found by other 
authors, both using a different form of cutoff ||^, ^ ^ ^ and those obtained 



using a sharp cutoff 0, ^ ||, |T2[. Ref.|T^ displays a table comparing results from 



different forms of the Local Potential Approximation. 

At 0{d^), foT N = 1 ■ ■ - 4, the results are already rather accurate for rj and 
z/, showing an improvement on the leading order results. This is also true of u 
ai N = 1, but for A^ = 2, 3, 4 the accuracy of uj is poorer than at (9°). For the 
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large cases however, uj is well estimated, v is not improved at 0(9^), while r] 
is dramatically underestimated - eventually by about a factor of 10. 

With the results exact at = oo, we see that there is a loss of accuracy in 
the approach to this limit. Indeed the approximation scheme appears 'too biased' 
towards the N = oo results. As we speculated in ref.0], this may be because all 
the appropriate fields have not been included in the effective action, in this regime. 
Indeed, it is known that at large A^, a massless bound state field also exists at 



the critical point p^ , |35[| . We should thus expect that a derivative expansion 
is not so well behaved, because the vertex functions are hiding within them the 
effects of this integrated out massless field. To ameliorate this behaviour, we 
should include the bound state explicitly as an effective 0{N) singlet field, then 
amongst the new set of fixed points in this enlarged space will be one with the 
same universal properties as the original A^ vector model, but with better behaved 
derivative expansion properties. We expect that similar considerations will apply 
to fixed points with fermions, particularly since the bound state fields here also 
correspond to the order parameter [a.k.a. fermion condensate) 
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A Numerical Methods 

In this section we outline the methods used to solve numerically our equations. 
We describe in particular, how to cast the equations manifestly as coupled second 
order differential equations, how the asymptotic expressions were used to set 
boundary conditions, and why it is necessary to develop these beyond leading 
order for K and V, and yet one further order for Z. We explain why relaxation 
was our method of choice for these problems, and why analytic forms for the 
integrals lij^k proved insufficient. We describe the singular behaviour seen as 
N ^ 1, and the relaxation from the exact results at A^ = oo that was used 
instead. We also point out some of the (simpler) issues that arose on solving for 
the perturbations. 

All the equations that we needed to solve numerically, form two point bound- 
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ary value problems. There are two methods for solving such equations: shoot- 



ing ^ and relaxation ||38[. Both methods were used, however relaxation is 
particularly suited to equations and boundary conditions involving complicated 
expressions which cannot be solved in closed form: there is no need to write the 
equation explicitly in the form dy/dx = f{x,y). It is also the best method when 
one needs to find solutions that depend upon some parameter, such as the value 
of A^. Once a solution for one value of N is found, one can use this solution as 
an initial guess for a close value of N. Given the nature of our problem it is not 
surprising that relaxation turned out to be our principal method. 

We formulate the equations as a set of non-linear, coupled second order dif- 
ferential equations, as follows [|19|. Noting that the equation for Z,cf. (p7D, [p^ . 



or for N = oo (pTl), contains powers of V", we differentiate the V equation 
( P^ ) and thus find an expression for V" in terms of V, V, V", K,K', Z and Z'. 
This is substituted back into the Z equation. Similarly, the (f'^ = boundary 
condition for the Z equation, contains powers of V"{0). We eliminated these 
in favour of V (0) , V (0) , K {0) and Z{0) by using the differentiated V equation. 
Similar transformations are required to turn the equations for the perturbations 
into manifestly second order differential equations. 

For the other set of boundary conditions we choose a value of y^^, which we call 
V^^gy, large enough such that the asymptotic expressions in appendix B become 
a good enough approximation to the solutions ||T^ (These expressions go 



beyond leading order. The leading order expressions (|38|) - pO|) are not sufficient 



as explained below.) We then use the asymptotic expressions for V,K and Z to 
provide boundary conditions for V, V, K, K', Z and Z'. The value of ipl^^ must 
not be chosen so large that numerical instability prevents us from obtaining a 
solution. There is a certain amount of trial and error in deciding where to set 
^asy'- cannot really decide where to set it until we know something about the 
solution. We checked that the solutions were stable against reasonable changes 
in the value of <^lsy 

It is necessary to develop the asymptotic expressions to beyond the leading 
order in order to find a solution []TU|. Thus for example in (PD|), if we just substitute 
the leading order results, viz. the first term in (^2]), for V^ifl^y) and V'{(pl^y), 
and attempt to solve for ^^"(v^asy) {^^ effect takes place in the numerical 
procedures), we see that there is no real solution for any finite v'asy Hence, we are 
forced to expand the asymptotic expressions to beyond leading order. A similar 
result holds for K^ip^). However, when we consider Z((y9^), we see that we need to 
take the asymptotic expressions to next-to-next-to-leading order: the correction 
to the leading order asymptotic behaviour does not involve Z" , hence the above 
procedure carried only to next-to-leading order leaves Z" incorrectly determined. 
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In order to ensure that Z" gets determined correctly (at least to leading order), 
we need to extend the asymptotic expression yet one further order. 

In the perturbed equations it is sufficient to use the leading order asymptotic 
expression |T^, As a check we calculated the corrections to the leading order 
behaviour and also used these. This made no difference to the values quoted. 



The solution of the equations is complicated by the fact that they are stiff |^ 



p8[ . This arises because small perturbations from the true solutions lead to a 
singular behaviour. Relaxation is particularly adept at handling stiff problems. 
In this way, the numerical solution of the leading order equations proved to be 
easy. It was only necessary to provide an initial guess and the relaxation routines 
quickly found the solution. 



The integrals (pq ) play an important part in the solution of the second order 
equations. It is possible to find a general analytic result for them, but these 
turned out not to be useful for the present purposes. The integrals are generally 
small numbers, but in the analytic expressions these small numbers often arise 
from the cancellation of large numbers. In this case roundoff errors play too large 
a part in determining their value using the analytic expression. Even performing 
the calculations using higher precision Fortran routines could not provide accu- 
racy beyond four significant figures in extreme (but realised) cases. Furthermore, 
whilst the equations (P^ ) - (^) contain about three hundred terms expressed in 
terms of the integrals li^j^k-, particularly for the larger j, k the algebraic expres- 
sions are so large that the equations expressed in closed form stretch to between 
twenty five and thirty thousand terms. The main effect of this dramatic increase 
is cumulative roundoff error. The attained level of accuracy proved totally in- 
sufficient and so numerical methods were used. We used an adaptive integrator 
from the NAG libraries. The problem was split into two parts: a ffist part from 
zero up to some finite value g^, and the remaining bit. To avoid roundoff errors 
we used an expression arising from asymptotic analysis to calculate the second 
part of the integral, and only used the numerical routines to compute the integral 
over the finite range ffist part. The point was set to be as small as possible (to 
avoid numeric round off in the numerical integration), whilst still allowing the 
asymptotic expression to be accurate. 

The numerical solution of the second order equation is then made difficult by 
the large number of integrals that need to be calculated numerically during each 
iteration. It was soon noted that the computational power required could not be 



provided by normal workstations. Parallel Fortran code, making use of MPI 
was developed to run on a sixteen node IBM SP2 system. This reduced run times 
for the relaxation code from up to one week to typically between ten minutes to 
one hour. 
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Initially we tried to solve the second order equations by relaxing away from 



= 1 using the numerical solutions obtained in [|T9| as initial guesses. However, 
this procedure failed because as — > 1 the fixed point solution becomes 
rapidly sharper at the origin, diverging as Z{ip'^) ~ in the limit = 1. 

This was confirmed by substituting = 1 and K = k — Zip^ [cf. (|4^) ] in the K 
equation (|35|), using the known results ||19[ for n and V , and solving for Z. This 
divergence is allowed because at A^ = 1 only the combination n = K + Zip^ has 
any meaning and this remains well behaved, as explained at the end of section 
4. Our successful strategy was to first solve the simpler N = oo equations, and 
then relax from 1/A^ = to finite A^ using these solutions as initial guesses. 

The A^ = oo, and K equations were solved analytically as described in 
section 5. The equation (^TJ) that determines Z at A^ = oo, was solved numeri- 
cally. Being first order, there is only one boundary condition, which is that the 
solution exists for all ip"^ and thus has the required asymptotic behaviour. For 
this problem we used the shooting method. An initial value of Z was used to 
integrate the equation out towards the second boundary using an eight point 
Runge-Kutta integrator. In general the initial value Z[0) would be wrong: if it 
was too large then the solution would tend to infinity at some finite y^^; if it was 
too small the solution would tend to minus infinity at some finite value of y?^. 
The correct value was then found by binary chop. It was particularly difficult to 
get past the point cp'^ = l/2v^. This is because at this point the right hand side 
of the equations has a 0/0 type behaviour, i.e. Z' = n^p"^, Z)/d{(p'^, Z), such that 
as (p"^ 1/2/ \/2, n{Lp'^,Z) — * and d(Lp'^,Z) — > 0, for the true solution. There 



are general methods of dealing with such singular points |3^, but our equations 
were sufficiently well behaved to avoid needing to use them. 

The shooting method was tried for 1/A^ > also, but this proved to be 
unsuitable. This is because to use shooting effectively we need to be able to reach 
the other boundary or have a particularly simple set of equations {e.g. the N = oo 
equations). The singularity structure of the equations makes it impossible to 
shoot from one side to the other and hence makes the more complicated problem 
unsuited to shooting. 

The main difficultly in solving the second order equation was removing the in- 
stabilities in Z near the origin. These manifested themselves in the form of sharp 
spikes in Z' near = 0. These were removed by writing Z{p>'^) = s{p'^)Z[ip'^), 
where s((/9^) is a known scaling function that makes as fiat as possible. 

This greatly speeded up the calculations and removed the spikes in the solutions. 
However even this could not quell the instabilities that arise as A^ — 1 where Z{0) 
actually diverges (as discussed above). As a consequence we found that accurate 
solutions even at A^ = 2 were much harder to obtain than those of A^ = 3 (say). 
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while the singular behaviour at = 1 prevented us from relaxing to = 0. This 
is the reason why we do not display = results at second order in table 1. 



In the equations for the perturbations we calculated the parts dependent only 
upon the fixed point solutions and stored them in a file, thus saving a large 



computational overhead |19, It was important that these functions, which 
are the functions that multiplied v, v', v", k etc. in each of the linearized equations, 
were well determined. The size of the z equation was such that small errors in 
the fixed point solutions made it difficult to determine these functions accurately, 
for large ip"^. This problem was solved by determining these functions using the 
asymptotic expressions for V, K and Z [cf. appendix B) and then matching them 
onto the ones calculated from the fixed point solutions. 

To conclude, we note that we checked that the eigenvalues were stable against 
reasonable changes in where ip^^^ was set, in the number of mesh points, and other 
numerical factors. A good indicator of numerical accuracy was to find numerically 
the redundant perturbation mentioned in section 4, and check that it matched 
the analytic form and that its eigenvalue vanished to good precision. (These 
properties are only completely recovered in the limit of an infinitely fine mesh.) 
We found we could numerically determine this perturbation and eigenvalue to a 
high degree of accuracy for all values of N. 



B Asymptotic expressions for K and Z 

From the second order equations (|3lD - (p7D, we determined the asymptotic be- 
haviour of V{ip'^) and K{ip'^) for large y?^ to next-to-leading order, and to 
next-to-next-to-leading order, for reasons explained in appendix A. With A^j, Ak 
and Az denoting constants, the results are 

( 1 V6(iV- 1) y/TT^ 1 x/6(l + r/) 

1,24 v^v^AT + 24 v^ATvO^TXyS^ 
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K{^^) - A, T^) + 2 (^v^ ( 1 + ) (4 r;M,3 + 16 r/^ A, 

+ 11 ry^ + 48 Aj^ + 55 .4^ 77^ A,"" - 52 ^3 

- 56 Ak" rf A, - 579 rj" A^ A,'' - 480 Ak^ rf A, - 80 77' 

- 144 rf" Aj^ - 703 77 Ak + 336 77 Ak^ + 40 77 

- 432 77 A,^ - 80 - 240 Aj" + 160 A^'' A) /( + A, 

{-AA, + 2A,7^ + A, + r^A,f{b-r^f'^NA,'l^) + 

^ 864 

^1 + 77 ((12 77^ AT + 28 77^ + 48 77 - 48 77^ - 32 77^ - 48 77) 
+ (6877^ ^ - 2477A^^ - 12 77^A^^ + 4477^^ + 877 ^ 
+ 12 77=^A^^ - 16 ^) Ak'^ + i^rf N A,'^ + ^r]N A,'^ + 21 rf A,'^ 
+ 69 77 ^2 + 66 rf A,'' + 24 A,'' + QrfN A,'') Au + brf Aj^ 
+ 15 77 ^ + 15^2 Aj' + bAj) /[A.'^'^^kN 

(( 277 - 4)^^ + 77^ + ^^)) ((^2)^-^) 



Z((^2) - yl^ ((^2-) I 1+^ J + /_ (^64 77^ AT - 65 77^ + 185 77^ - 172 77^ AT - 716 77 

+ 652 rf + 604 rf N - 452 rf - 2524 7; + 2348 r]N - 1184 

+ 1264 A^) /(^yXiV (-2 + 77)2(l+77)(5-r7 )^/^) 
^ x/6 (2077^ - 33 7/3 - 67/^ - 587/ - 24) ( AT - 1 ) X 

24 (l + ;^)3/2(_2 + ^)2JV^ ' 
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